Effect of dipolar moments in domain sizes of lipid bilayers and monolayers 
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Lipid domains are found in systems such as multi-component bilayer membranes and single com- 
ponent monolayers at the air-water interface. It was shown by Andelman et al. (Comptes Rendus 
301, 675 (1985)) and McConnell et al. (Phys. Chem. 91, 6417 (1987)) that in monolayers, the 
size of the domains results from balancing the line tension, which favors the formation of a large 
single circular domain, against the electrostatic cost of assembling the dipolar moments of the lipids. 
In this paper, we present an exact analytical expression for the electric potential, ion distribution 
and electrostatic free energy for different problems consisting of three different slabs with different 
dielectric constants and Debye lengths, with a circular homogeneous dipolar density in the middle 
slab. From these solutions, we extend the calculation of domain sizes for monolayers to include the 
effects of finite ionic strength, dielectric discontinuities (or image charges) and the polarizability of 
the dipoles and further generalize the calculations to account for domains in lipid bilayers. In mono- 
layers, the size of the domains is dependent on the different dielectric constants but independent 
of ionic strength. In asymmetric bilayers, where the inner and outer leaflets have different dipolar 
densities, domains show a strong size dependence with ionic strength, with molecular-sized domains 
that grow to macroscopic phase separation with increasing ionic strength. We discuss the impli- 
cations of the results for experiments and briefly consider their relation to other two dimensional 
systems such as Wigner crystals or heteroepitaxial growth. 

PACS numbers: 82.45.Mp,87.16.Dg,87.14.Cc 



I. INTRODUCTION 

Single component monolayers of phospholipids or fatty 
acids at the air-water interface exhibit stable, micron 
sized domains of coexisting phases [J H, E|. Such do- 
mains have not been observed in single component bilay- 
ers, but vesicles consisting of mixtures, such as ternary 
mixtures of phospholipids, sphingolipids and cholesterol, 
which have been the subject of intense investigations 
[3, Q, do exhibit stable domains. Most of the natural 
occurring phospholipids and sphingolipids are zwitteri- 
onic @ and exhibit a permanent electric dipole moment, 
which is quite significant, around 18D[7|. In a monolayer, 
these dipoles all point in the same direction, with its ma- 
jor component perpendicular to the air-water interface. 
In single component bilayers, the permanent dipoles in 
both leaflets point in opposite direction and the mem- 
brane does not have a permanent dipole. In many im- 
portant situations, however, the lipid composition within 
the inner and outer leaflet is asymmetric. In most cell 
membranes for example, phosphatidylethanolamine (PE) 
and phosphatidylserine (PS) are preferably found in the 
inner leaflet (in contact with the cytoplasm) while sph- 
ingomyelin (SM) and phosphatidylcholine (PC) are more 
abundant in the outer leaflet § leading to a permanent 
dipole density for the membrane, similarly as in mono- 
layers. In fact, it is well established that a difference in 
electric potential across a membrane due to dipolar po- 
tentials is quite significant, between 100 mV and 400 mV, 
and that this potential difference is extremely important 
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for a variety of biological processes [9( . 

What is the effect of domains contained in a matrix 
with a different dipole density? In a series of papers 
Andelman et al. flO . Illll and McConnell and collabo- 
rators [H, EH EH Il5l . Il6j | have provided theoretical and 
experimental evidence that dipolar domains in monolay- 
ers have a profound effect in the phase diagram, as the 
electrostatic free energy cost of assembling such domains 
limits its size, preventing the monolayer from reaching 
complete phase separation. In their original calculation, 
McConnell and collaborators did not consider several ef- 
fects that are potentially relevant, such as a finite ionic 
strength, the inclusion of image charges and the polariz- 
ability of the dipoles. Those effects were considered by 
Andelman et al. [ll| for geometries consisting of stripes. 
Furthermore, in view of the importance of lipid domains 
in membranes, it is of great interest to analyze how dipole 
density inhomogeneities may affect the phase diagram 
[ijj, similarly as it is the case for simple monolayers at 
the air-water interface. In this paper, we generalize the 
calculation by McConnell and collaborators for circular 
domains to include finite ionic strength, image charges 
and the polarizability of the dipoles and further extend 
the results to dipolar domains in bilayers. 

The phase diagram of lipid mixtures is a result of the 
microscopic characteristics of the lipids involved such as 
hydrocarbon chain and length as well as degree of satura- 
tion, possibility of forming hydrogen bonds, etc.. It is ex- 
pected that the effect of these interactions is parameter- 
ized on longer scales by quantities such as a line tension 7 
between coexisting domains, a spontaneous curvature cq 
of the different phospholipid species, the membrane sur- 
face tension r or its bending rigidity k (see for example 
[l8l . [l9j and references therein) . The incorporation of the 
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FIG. 1: An asymmetric bilayer with a domain with a dipole 
density. 



microscopic permanent dipoles of the lipids in these long 
wavelength descriptions is subtle. As an example, let us 
consider a bilayer containing a circular domain with non- 
zero dipolar density within a matrix of an hypothetical 
dipole free lipid, as shown in Fig.[TJ The electrostatic free 
energy cost of such domain can be estimated in analogy 
with the monolayer case [12|, ll3j 



F « R(a + blogR) 



(1) 



where R is the size of the domain. If the coefficient b 
were vanishingly small, the net electrostatic cost of the 
dipoles grows with the perimeter of the domain and may 
be lumped into a "renormalized" line tension contribu- 
tion. For the case shown in Fig. [IJ however, the coeffi- 
cient b is not negligible, and the free energy contribution 
of the permanent dipoles must be considered explicitly. 
In some other situations, however, the coefficient b may 
vanish. It seems plausible that if the domain in Fig. [T] 
was opposed by an identical domain of the same lipid, 
the membrane would not exhibit a permanent dipole mo- 
ment and the electrostatic contribution would consist of 
a short-ranged interaction that might be incorporated 
into the line tension. Furthermore, ionic strength, by 
weakening the electrostatic interactions, can also influ- 
ence domain sizes. In this paper we present a detailed 
quantitative analysis for these questions. 

In all previous treatments of dipolar domains (such 



as 



131 ] and all subsequent work), the free energy has 
been computed by directly summing over dipoles. In 
this paper, we will consider two oppositely charged sheets 
and obtain the free energy of the dipoles as the first 
non-trivial coefficient in the limit of small sheet sepa- 
ration. Obtaining the dipolar result in this way is partic- 
ularly convenient for implementing the boundary condi- 
tions that should be met by dipole densities within layers 
with different dielectric constants. 

The organization of the paper is as follows. In Sect. |TT] 
we derive an explicit form for the free energy of the 
system. A summary of the exact solution for the elec- 
trostatic free energy of dipolar domains is provided in 
Sect. IIII1 The details of the solution are mathematically 
dense so we defer the complete derivation to appendix [A"l 
In Sect. IIVI we analyze the theoretical implications for 
monolayers and bilayers domains. We conclude our pa- 



per with a discussion on some of the consequences of our 
results in Sect. El 



II. ELECTROSTATIC FREE ENERGY 

We consider a system of three slabs with different di- 
electric constants and Debye lengths. Two disks of op- 
posite surface charges parallel to each other are located 
on the boundaries of the middle slab, as shown in Fig. [5] 
The distance D between the two oppositely charged disks 
is assumed smaller than any other length in the problem 
and therefore the two disks represent a dipolar density. 
In the dilute electrolyte regime, the free energy of the 
system consists of an ideal and an electrostatic term, [2(| 

F = k B Tj2 J d 3 rn a (r)(log{vn a {r)) - 1) + 



— / d 3 rE ■ D 
8tt 



(2) 



where n a (r) = exp(— 9a ^^ ) is the number density 
of ions of type a, ionic valence q a and bulk concentra- 
tion . It is convenient to express the free energy in a 
different form. For the ideal part we write 

/ d 3 ™ Q (r)(log( vn a (r)) - 1) 

a J 

= "E^T / rfWr)<&(r) + 5> a (log(nf«)-l) 

a ^ a 

= - E U / d 3 rn a (r)<f(r) + £ N a Qog(v^) - 1) - 



5>f j d 3 r(exp(- 



gqe$(r) , 
k R T ' 



(3) 



where the last term follows from the identity N a = 
J d 3 rn^ exp(— q "^^ ) , where N a is the number of ions 
of type a within the solution. The electrostatic term is 
written as 



— f d 3 rE D= — - [ d 3 rV$(r) • D (4) 
8n J 8ir J 

i J d 3 re$(r)^g Q n a (r) + i J d 3 re$(r)a(r) 



The previous expressions are general for an interface in 
contact with a dilute ionic solution. Dipolar interactions 
are weak so we further assume the "Debye-Huckel" ap- 
proximation Tpjf? << 1. The ion number density be- 
comes 

n a (r) =„„(!- & _ + - {kBT)2 + ..). (5) 
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Upon using the charge neutrality condition q a n^ = 0, 
it follows 

£ «.«.(') = +•■•(«) 

Inserting Eq.[3J Eq.H Eq.EJand Eq.[5]mto Eq.Hthe free 
energy is 

F=±j d 3 r<P(r)a(r)+k B Tj2N a (H^v)-l). (7) 

The first term is the electrostatic free energy of the 
charged disks. It is remarkable that both the entropic 
and the electrostatic terms involving the ions explicitly 
cancel out. The second term is the ideal free energy of 
the ion species in solution. The subsequent analysis will 
be performed within the canonical ensemble where the 
total number of particles of the different species is fixed, 
so the ideal free energy is constant and will be ignored 
herein. For uniformly charged interfaces, the free energy 



can be further simplified to 

f°° D D 

F = ir\a\j dpp(<P(p,--)-<P( P: -)). (8) 

A complete evaluation of the free energy requires the de- 
termination of the electric potential, which follows from 
the Poisson-Boltzmann (PB) equation 

V 2 $(r) = -L$(r) (9) 

where is the Debye length. In regions where no ions 
are present (Xd — oo) and the previous equation reduces 
to the standard Poisson equation. In this paper, we con- 
sider problems with cylindrical symmetry. The solution 
to the PB equation is 

$(p, z) = a{k) exp(±fcz) J Q {p{k 2 - l/X 2 D ) 1/2 ) , (10) 

where Jq{x) is the 0th order Bessel function and k indexes 
different solutions. 

The electric potential is 



/ °° dkJ (kp) 

$(j>,z) = { rdkj (kp)( 



cxp(-(fc 2 + l/(A^) 2 ) 1/2 *) 

(fc2 + l/(AJ,)2)l/2 

cxp((fc 2 + l/(A b n ) 2 ) 1/2 *) 



a 2 (k) - 



(fe 2 + l/(Ap 2 )i/ 2 

/o dkJo(kp) P ( fc2+1 /( A c^2)i/2 a,i(k) 



z > D/2 

^(^/(©iy^ asW) -D/2 <z<D/2 

z < -D/2 



(11) 



The condition that the electric field must vanish at in- 
finity (E z (oo) = E z (—oo) = 0) has been already imple- 
mented. The values of A^ b,c are the Debye lengths in each 
region. We further assume a different dielectric constant 
£ a,b,c f or ggyjjj re gi on . 

The coefficients {a,i(k)}i—i...4 are determined from the 
boundary conditions 

Hp, (d/2)+) 

$(p,(-£>/2)+) 
e a E z (p, {D/2)+)-e b E z { Pl {D/2)-) 
e b E z {p, (-D/2)+) - e c E z {p, {-D/2)-) 



where E z is the component of the electric field along the 
z-axis, which is assumed to be perpendicular to the plane 
defined by the disks. The boundary conditions defined by 
Eq. rT3] are easily implemented by recalling the following 
identity for Bessel functions 

@{R-p)= dkJo{kp)J!(kR) , (13) 
Jo 



where O(x) is Heaviside-0 function defined by O(x) = 1 
if x > and <d{x) = for x < 0. 

For further reference, we quote the free energy of two 
infinite disks of uniform and opposite charges filled with 
a medium of dielectric constant ej. The free energy of 
this system is equivalent to the elementary formula of an 
infinite capacitor 



F GC = ^l|a| 2 i? 2 A D . (15) 

We note that upon replacing ej by the dielectric constant 
of the solvent e w and the finite size D by the Debye length 
\d , the "capacitor" formula Eq. rjj] is identical with the 
previous equation. With the explicit form of the potential 
Eq.[TT]introduced into Eq.0 the free energy of the system 



= <£{p,{D/2)-) F c = —\a\ 2 R 2 D. (14) 

= <f{p,{-D/2)-) 

= — 47r|<7|6(i? — p) We also will need the non-ideal free energy for the 
= 47r|cr|0(i? — p) Guoy-Chapman diffuse layer of a weakly charged system 
Ac/Ag << 1, where A^ is the Debye length and Xq is 
^ ' the Guoy-Chapman length. The free energy is given by 
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£j -ft Ad 



(16) 



where T T {x, y) is a function specific for each situation, 
which according to Eq. [TJ] satisfies ^ rX (0, 0) = 1. In the 
following we provide an explicit calculation for the func- 
tion ^{x^y) in different cases. 



III. EXACT RESULTS FOR THE FREE 
ENERGY IN SEVERAL CASES 

We provide an explicit solution for the following three 
problems shown in Fig. [2] In all cases, two infinitely 
thin disks of finite radius R and uniform but opposite 
charges are separated by a distance D, and it is assumed 
that R » D. In Problem A) the two disks are on the 
boundary of a medium of dielectric constant £/ , with only 
one domain in contact with a salty solution with solvent 
of dielectric constant e w (water) and the other domain in 
contact with a medium of dielectric £a- In problem B) 
both domains are in contact with a salty solution with 
solvent of dielectric e wi while the interior contains a po- 
larizable medium of dielectric constant £/. In Problem 
C) the two domains are submerged in a solution with 
solvent of dielectric e w . 

The free energy of these three systems can be com- 
puted exactly. The most relevant results are summa- 
rized below. Explicit formulas and complete details of 
the derivation are provided in appendix lAl 



A. Finite circular dipole domain in a medium of 
dielectric ei with one side in contact with a salty 
solution and the other with a medium of dielectric 

£A 

The free energy for this system is given by the general 
expression 



271" , ,9 9 , 

= \a 2 R 2 D(l- 

si 

tth ad ad 



(17) 



where /i(y) and f^iy) are analytic functions of y. In 
the limit y — > (R/Xd << 1, low screening limit) and 
y — > oo (R/Xd » 1, high screening limit) both f\ and 
f% can be computed analytically. The results for /i are 



/i(0) = 2 

AM = 2 



E w £a 



ei(e w + e A ) 
£i 



(18) 



Explicit formulas for f% and further details of the calcu- 
lation are provided in appendix IA 11 
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FIG. 2: Two domains of uniform and opposite charges in 
contact with: A) one side with a salty solution, a medium in 
between and a different medium in the outside. B) Both sides 
in contact with a salty solution and a medium in between. C) 
submerged in a salty solution. 



B. Finite circular dipole domain within a medium 
of dielectric ei both sides in contact with an aqueous 
solution 

The free energy adopts different forms. In the low 
screening limit it reduces to a similar expression as in 
the previous case, 

F b = ^l^R^ _ JL {fl HR/D) + h)] (19) 
£/ nR 

where /i = — and the explicit expression for /2 is quoted 

in Eq. rsa ei 

In the high screening limit, there are two situations 
that need to be considered, R » Xp >> D and R » 
D » X D . For R » X D » D it is found 

F B = 2—\a\ 2 R 2 D(l + 0(D/X D )). (20) 
£/ 

This result is derived for a simplified problem discussed 
below. In the limit R >> Xd >> D it follows 



,2tt 2 



F b = 2— \a\ 2 R 2 X 



D ■ 



(21) 



Explicit formulas and detailed calculations are provided 
in appendix IA 21 



5 




Coarse Graining 




(Uniform density of charge) 



Hydrocarbon 
/ Tail (e HC ) 



^ Head group 
^region 

Water (e w ) 



FIG. 3: A Langmuir monolayer with a domain with a finite 
dipole density. The monolayer head group is represented as a 
medium with dielectric constant ej within two opposite charge 
planes with uniform surface charges, one side in contact with 
an aqueous solution and the other with a medium of dielectric 
constant shc (the hydrocarbon tail). 



C. Two oppositely charged disks in solution 

The low screening limit has the same form as in Prob- 
lem A. In the high screening limit, and for D >> \r>, 
the system consists of two infinitely thin disks with dif- 
fuse layers on both sides interacting through a screened 
Coulomb potential, 



2n 2 



a\ 2 R 2 \ D (l-cx V (-D/\ D )). (22) 



In the high screening limit, but for R >> \n >> D, the 
leading term of the free energy is the "capacitor" term 
and is given by 



— \a\ 2 R 2 D{\-—) 

Ej ZAd 



(23) 



The first correction is independent of the radius R of the 
domain. The details of the calculation are provided in 
appendix IA 31 



coarse-graining suggested in Fig. [3l which consists of ap- 
proximating the head group region as a medium of dielec- 
tric constant ei and the hydrocarbon tail as an infinite 
medium of dielectric constant Ehc- From the results in 
Subsect. lTTTAl it is 



F tot = (27r 7 i? + F J 



ttR 2 



(24) 



where 7 is the line tension between the two phases. The 
free energy explicitly depends on the domain radius R, 
and is minimized under the constraints of both constant 
dipolar density /1 and minority phase area A s . The equi- 
librium radius of the domains R eq is given from 



R eq = Dexp(l - ^)exp(-^n) ; 



£exp(i% (25) 



The above expression is only exact in the limits y — > 
(low screening) and y — > 00 (high screening), as otherwise 
terms involving derivatives of f(y) are present. From 
Eq. [TH1 as a result of e w >> ea = shc, it is found 
/i(0) « /i(oo). It also follows from Eq. IS] that / 2 (0) « 
/2(oo). Therefore, for both the low and high screening 
limits, the domain radius is the same and is given by 



R eq w £>exp( 



2 

£/7 



-)• 



(26) 



The term in the exponent is only half of the value origi- 
nally obtained by McConnell [12, HH and it is multiplied 
by the factor e 2 /ehc- For the case e = 1, it agrees, 
upto the 1/2 factor, with the result for stripes obtained 
in [ll[ . The origin of the 1 /2 factor is the discontinuity in 
the dielectric constant between air and water, while the 
presence of the different dielectric constants is due to the 
finite polarizability of the dipoles and the surrounding 
medium. 

The same formula applies for a more general situa- 
tion where the matrix has a non-zero dipole density [i\ 
and the domains a dipolar density fi2- In that case, the 
quantity appearing in the exponent is the dipolar density 
difference \i = — 



B. Bilayers 



IV. IMPLICATIONS FOR MONOLAYERS AND 
LIPID DOMAINS 

A. Monolayers 

We consider a monolayer with dielectric constant ei in 
between water e w w 80 and air ea%t ~ 1- The monolayer 
consists of a minority phase of n s domains with radius R 
and dipolar density fi = aD contained within a matrix 
of a majority phase with zero dipolar density, as shown 
in Fig. [3J The area spanned by the minority phase is 
A s = n s irR 2 . The free energy is constructed from the 



We now consider a bilayer with a circular domain, 
where the dipolar density of the outer leaflet (/ii) and 
the inner leaflet {1x2) are different. Similarly as in the 
case for monolayers, the dipolar densities are to be in- 
terpreted as the excess dipolar density from the matrix 
reference state. The domain has radius R and the verti- 
cal separation between the two leaflets is L, as shown in 
fig- HI The bilayer is surrounded by water (e w ps 80). As 
shown in appendix [Bl the electrostatic free energy cost 
is of the form 

F e i e - -(mi - fi 2 ) 2 Rln& + "f'2irR, (27) 
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FIG. 4: Dipolar domains in a bilayer. In a first approxi- 
mation, the two leaflets of the bilayer are modelled as two 
monolayers, an inner leaflet monolayer with dipolar density 
fi2 = 02-D and an outer leaflet monolayer with /ii = aiD, 
following the same steps as in Fig. [3] At longer distances, 
the net dipolar moment of the two monolayers can be repre- 
sented as if were originated by two oppositely charged disks 
with charge density a' — (cti — (J2)D/L, with an effective di- 
electric constant ei taking into account both hydrocarbon tail 
and head groups polarizabilities. 



this configuration leads to the same dipolar moment. It is 
expected that higher order moments provide sub-leading 
corrections, which nevertheless may be incorporated as 
any other short-range interactions, if necessary. The di- 
electric constant ei is an effective value that includes con- 
tributions both from the hydrocarbon tails and the head 
group and therefore it should be expected that ei ~ 4. 
The different steps in the coarse graining process are illus- 
trated in Fig. |U Under these assumptions, the problem 
of an asymmetric domain is reduced to the problem dis- 
cussed in Subsect. IIIIBI In the low screening limit, the 
radius of n s domains is computed as for the monolayer 
case, and is given by 



Re 



L exp( 



e?7 



i)> 



(28) 



where e / is the effective dielectric constant including con- 
tributions from the hydrocarbons and the polarizabilities 
of the dipolcs. This formula differs from the analogous 
result for monolayers in the factor s w /2ehc ~ 10 in the 
exponent, which implies that if a given dipole density and 
line tension in a monolayer leads to an equilibrium radius 
R eq ~ 10/iM , the same dipolar density and line tension 
in a bilayer results in much smaller domains, of radius 
comparable to the molecular scale (of the order of L). 
Alternatively, the previous formula shows that the con- 
dition for large (micron sized) domains to exist is given 
by 



£ w fJ- 2 



8 -> Ap « eiA 



8e w L 



(29) 



where terms proportional to R 2 have been ignored as they 
are not relevant for this argument. For a symmetric do- 
main (/ii =1-12) there is no logarithmic contribution and 
the electrostatic free energy is of the same form as the 
line tension. The above result has been derived without 
including the effects of ionic strength or different dielec- 
tric constants, but we assume it true under these more 
general conditions. 

If the membrane contains an intrinsic asymmetry be- 
tween the inner and outer leaflets (for definiteness, we 
assume [i\ > [12), there is a net dipolar moment. We 
consider domains much larger than the vertical separa- 
tion of the two leaflets (R >> L), and construct the free 
energy by first coarse-graining each separate leaflet of the 
bilayer following the same steps as in the monolayer case. 
The bilayer has a net dipolar moment, which we model 
as if it were produced by two disks of uniform but oppo- 
site surface charge density a', where a' = (ai —o<i)DnlL 
and Dh is the width of the head group region. This ar- 
gument should not be understood as implying that the 
coarse-graining of an asymmetric bilayer should lead to 
equal and opposite charge distributions, but rather that 



where p is the permanent moment of the phospholipid, 
Aq is the molecular area and we assumed that the line 
tension is 7 « *g£. With A « 60 A 2 , V « 10 A, 
the formula yields Ap << O.lej (D). If we assume £/ ss 
4 then Ap << 0.4 D. That is, a net dipolar moment 
difference of the order of 0.4 D or less is needed for large 
(micron sized) dipolar domains to be observed. Dipolar 
moments for phospholipids are much larger, of the order 
of 15 D 7], which implies that domains with significant 
asymmetry will be sub-micron in size. 

In the high screening limit, the results are qualitatively 
different. When the Debye length becomes of the order 
of the domain radius, the electrostatic interactions be- 
come short-ranged and the free energy, which is given 
in Eq. 12 1 1 does not show the logarithmic contribution 
present in Eq. \%7\ In that case, the thermodynamic sta- 
ble state of the n s domains is completely dominated by 
line tension and the free energy is minimized by a single 
circular domain (assuming that the resulting line tension 
is positive) of the minority phase within a matrix of the 
majority phase. A summary of the different results is 
shown in Fig. [5] 
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FIG. 5: Schematic representation of the different systems 
studied upon variations in ionic strength. 

V. DISCUSSION AND CONCLUSIONS 

We presented a detailed calculation showing how the 
competition between line tension and electrostatic dipo- 
lar interactions determines the size of lipid domains, thus 
extending previous results by McConnell and collabora- 
tors fl2l . Il3| by including finite ion strength, variations 
in the dielectric constant and the polarizability of the 
dipoles and further extending the results to account for 
dipolar domains in lipid membranes. The main conclu- 
sions are summarized in Fig. [5j for lipid monolayers at 
the air-water interface it is found that the size of the 
lipid domains is independent of the ionic strength, but 
is sensitive to the discontinuities in dielectric constant 
and the polarizability of the dipoles. In lipid vesicles, 
it has been shown that for symmetric domains, where 
the net dipole density in the outer and inner leaflet is the 
same, the electrostatic contribution to the free energy can 
be absorbed into the line tension coefficient. For asym- 
metric domains, which have a net dipolar density, it has 
been shown that the large dielectric constant of water 
is responsible for a dramatic decrease of the size of the 
domains as compared with the size obtained for a mono- 
layers with the same line tension and dipole density. It 
has also been shown that small variations in dipolar den- 
sity (2% or more from the equilibrium dipole density of 
a zwitterionic phospholipid) are enough to limit domain 
sizes to sub-micron scales. Opposite to the situation with 
monolayers, dipolar domains in bilayers are sensitive to 
solution ionic strength. At small Debye lengths, of the 
order of the domain radius at zero ionic strength, dipo- 
lar electrostatic effects gradually become negligible and 
the system is dominated by line tension, resulting in a 
single large circular minority domain if the line tension 
is positive. Our study has been focused on the limiting 
cases of high and low ionic strength, but exact analytic 
expressions have been provided for the free energy and 



the electric potential. Problems that have not been dis- 
cussed in this paper, such as values for domain radius at 
intermediate ionic strength, ionic distributions or poten- 
tial differences across the interface also follow from these 
expressions. 

Our results show the subtle role played by the intrinsic 
dipoles in lipids membranes. At high ionic strength, the 
effect of dipoles is short-ranged, but at low ionic strength 
dipolar interactions are only negligible for symmetric do- 
mains, and in any other situation, dipolar interactions 
must be included explicitly, either in continuum models 
0, Ell or in numerical simulations [12] • 

Recent experimental and theoretical studies have stud- 
ied the interactions of circular and planar dipolar do- 
mains [23[ and interactions within two domains [24j. As 
compared with previous theoretical studies [ItJ of domain 
shape transitions in lipid bilayers, our study includes the 
effects of a net dipolar moment and the dielectric con- 
stant of water. More importantly, it shows that at high 
ionic strength the effect of dipolar interactions is short- 
ranged. Other theoretical work [25[ has discussed in- 
terfaces within the framework of the Debye- Huckel ap- 
proximation, the same approximation used in this study, 
but has been focused on effective interactions between 
charges, so these results are somewhat complementary 
to the ones presented here. 

A number of situations have not been considered in 
this paper. Our study has been restricted to circular do- 
mains, while other shapes, such as stripes, are also of in- 
terest [T3, [2^] . We further assumed a planar membrane, 
thus neglecting the effects of Gaussian curvature, which 
has been shown that, at least for solid domains, may lead 
to significant effects [H, [2?], HI] . The results presented 
have ignored domain-domain interaction and i ts p ossible 
superstructures. This problem is discussed in [291 ]. where 
it is shown that triangular lattices of circular domains are 
usually found, and a discussion is presented on the tran- 
sition from circular to other domains and other phases 
such as stripes [301 ] . 

Problems where dipolar-like domains, similar to the 
ones considered in this paper, have been discussed in 
other contexts. The free energy describing the phase 
transition between a Fermi liquid and an insulating 
Wigner crystal phase in a clean two dimensional elec- 
tron can be mapped to a dipolar interaction and a line 
tension term. After a detailed analysis, the authors in 
[3l| showed that a direct first order phase transition can- 
not occur, and instead, it was shown that the transi- 
tion should proceed as a sequence of intermediate "mi- 
croemulsion" phases, consisting of domains of one phase 
within a matrix of the other phase. Our results are in 
complete agreement with [3l| , but show that a direct first 
order transition can be obtained as a function of the De- 
bye length, if the electrostatic interactions are screened. 
Previous MD simulations [H], [33[ have shown that transi- 
tions from triangular superlattices to striped domains are 
a general feature resulting from increasing the strength 
of a short-ranged attraction force with respect to a long 



8 



range repulsion. Upon identifying the short-range attrac- 
tion with the line tension and the long range repulsion 
with the dipolar interactions, these results are qualita- 
tively similar to the problems studied in this paper. It 
seems therefore that "microemulsion" phases are a com- 
mon feature of systems near a phase transition consisting 
of a long range repulsion and short a range attraction. 
In fact, this has been shown rigorously for general 2D 
Coulomb systems |35j |. Another interesting example is 
provided by heteroepitaxial growth of 2D islands in solid 
surfaces, where the elastic free energy cost of assembling 
an island also contains a logarithmic term (completely 
analogous with the dipolar case studied in this paper) 
and competes against the line tension. In [36(, it was 
shown that the addition of domain-domain interactions 
in the form of dipolar interactions leads to the formation 
of triangular superlattices. These results show that in- 
teractions among finite circular domains are responsible 
for the formation of superlattice structures of circular do- 
mains, whose size can be computed from balancing elastic 
(or dipolar) free energy against line tension, and provide 
yet another example illustrating the expected general re- 
sults of phase separations in systems with long range in- 
teractions. 

It is quite interesting to compare the exact solutions 
discussed in Sect. IIIII with the corrections to a finite size 
capacitor, as discussed Landau & Lifshitz [37]. Al- 
though the charge distribution in the plates of a finite 
size capacitor is not uniform, as we assume in our calcu- 
lations, the coefficient of the logarithmic term fi in the 
electrostatic energy for sj — ea — £w, Eq. IA7I is iden- 
tical with the analogous formula quoted in Landau & 
Lifshitz H3. 



This paper has shown the subtle role played by dipolar 
moments. We showed that the zwitterionic nature of 
the phospholipids leads to new effects both in bilayers 
and monolayers. It is our expectation that the results 
presented will be useful in other areas. 
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APPENDIX A: DETAILED DERIVATION OF 
THE EXPRESSIONS FOR THE FREE ENERGY 

1. Finite circular dipole domain in a medium of 
dielectric ei with one side in contact with a salty 
solution and the other with a medium of dielectric 

£A 

This situation corresponds to case A) in Fig. O There 
are three dielectric constants, the medium in between 
the disks (e = £/), the medium outside the disks (e c = 
£a) and the dielectric of the solvent e° = e w (water). 
The electric potential in the different region follows from 
Eq. [TO] with both \ b D = \ c D = oo. 

The coefficients {a,i(k)}i = i.± are determined from the 
boundary conditions Eq. [12] with the help of Eq.[l3l The 
explicit result for the 04, 03 coefficients are 



-(sinh(fc£>/2) + §f cosh(kD/2)) S inh{kD/2) exp((k 2 + 1/ X^) 1 / 2 D /2) 8 ir\a\R „ , 

a>i\k) — , ,, „ ; Ji(kR) 

(fca+1 '/ AD)1/a )+smh(fcJ/2)(^ + gp^faiT* ) ^ 

( ( gB + 77 ™sh(kD/2)) sinh(fc£/2) exp((fc 2 + l/X%)^D/2) ^MR 

aAk) — ... " , Ji(kR) 

y ' cosh(fc£>/2) , E ,„ ea \ 1 ~j n hf/ r 1 1 1 kei 

ei \ k + (fe2 + i/ A£ ,)i/2 ) + Sma(KU/Z)(-^j- -t- 2 + 1/x } i/2) 1 



r 



The explicit form for free energy follows from Eq. [8] where T A {x, y) is given by 

F = \a\ 2 R 2 Df A (-, — ) , (A2) 

£j it ad 



(Al) 



A(x . = f°° dk j2 sinh(fcx/2) cosh(fc a: /2)(f + (k2+y2)1/2 ) + sinhjkx) {k2+y2)1/2 

Jo k 2 lU x cosh(fc a: /2)(| + {k2 + b y2)1/2 ) +smh(kx)(f + (fc2+ * 2)1/2 ) ' 

I 

In order to alleviate the notation, we have introduced new variables a — £ w / £1 and b = Sa/sj. In the strict R » D 
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limit, the free energy corresponds to two infinite circular 
plates, as the function T A satisfies 



lim^(x,2/) = l 



(A4) 



irrespectively of the value for y. The first correction 
0(D/R) to the free energy is of the form 



F = — \a\ 2 R 2 D(l 



^{ h (^-)\u{R/D) + h(^-)}) 



where fi(y) and / 2 (2/) are two analytic functions, whose 
limiting expressions in the low and high screening limits is 
provided below. The general derivation of the expression 
Eq. lA5l follows from the same steps as the formula derived 
below for the restricted case of the low screening limit. 



a. The low screening limit (R « \d ) 
The function T A must be considered in the limit y 



(A5) 0, 



T A (x,0) 



4 f 00 dk 
x J Q k 2 

2 f°° dk 



x 



Jl{k)sin\i{kx/2) 



(a + b) cosh(fcr/2) + 2 sinh(fcr) 
(a + b) cosh(fca;) + (ab + 1) sinh(fca;)) 



o 



— Jf(k)(l-ex P (-kx))- 
k z 



(2-a-b) 
2(a+b) 



(1 — exp (— kx)) 



(ab+l-a-b) 
2(a+b) 



(1 - exp(-2fcx)) 



r 



(A6) 



This function is of the general form Eq. IC1[ where A 



2-a-b 



and B 



lo+l — a— b 



The general expression for the 



2{a+b) 2(a+b) 

asymptotic behavior follows from the result in Eq. [C2j 



/i(0) 
/a(0) 



= 2 



ab 
> 

J a + b 
ab 



£ W £A 



a + b 
2(a + b) 



HS/(e 1/2 )) 



e A ) 



(2 - a - b)(2a + 2b) 



ab+l- 
4 -2a 



ab+l 



(a + 6) 2 

- a — b , 



2b _ .ab+l 

— Q2{ ^bTT 



ab+l + a + b 

a — 6 N 



a 



(A7) 



where the functions Gi are defined in Eq. lC3l If both salt- 
free mediums have the same dielectric constant sj = £a, 
then b = 1 and the expression for fa simplifies to 



/ 2 (0) = HS/e 1 / 2 ) - 



1 



2(1 



ln(4/e). (A8) 



In [2lJ the coefficient /i(0) for a dipolc domain with 
no polarization si = 1 was quoted, and the result is in 
agreement with Eq. IA7I McConnell and collaborators 
[l2l Il3j have computed the free energy in the limit e w = 
ea — — 1 by integrating the individual contributions 
of single dipoles and introducing a short-distance cut-off 
5 to avoid a self-energy singularity. Their result is 



A(0) = 1 

/ 2 (0) = He 2 S/(AD)) , 



(A9) 



which is identical with Eq. IA7I if the cut-off is taken as 
S = D/{2e 3 ^ 2 ). 

b. The high screening limit R >> Ad 

We now consider the limit where the Debye screening 
is much shorter than the radius of the dipolar domain. 
We therefore need to consider the limiting case y — > 00 
in the scaling function Eq. IA3I The resulting T A is 



2 r dk TV^h (uk 1-1(1 -exp (-fee)) 

_ / I2 ^i W(l-exp {-kx) - 

X In k Z 



l + ^(l-exp(-2fcr)) 



(A10) 



which is of the general form Eq. IC1[ with A = —1/2 and coefficients follow as 
B = (b — l)/2. The expressions for /i(oo) and /a(oo) 



/i(oo) = 26-2^ 
si 

/ 2 (oo) = 2&ln(8/e 3 / 2 ) 



+ 



(All) 



,6 — 1, 
b + l\ g ^b+-i )+ 
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Both /i(oo) and /2(oo) are independent of e w . In the 
high screening limit, the electric field does not penetrate 
in the water, and therefore no reference to the dielectric 
properties of the water should be expected. 



domain. The potential is screened everywhere except for 
the thin layer in between the two disks, where the dielec- 
tric constant is ej. The electric potential follows from 
Eq. El with X a D = X° D = X D and X b D = oo. 



2. Finite circular dipole domain within a medium 
of dielectric ej both sides in contact with an aqueous 
solution 

Upon imposing the boundary conditions Eq. [12] the 
This situation differs from the previous case in that coefficients satisfy ai(k) = —a^k) and a,2(k) = —az(k), 
the aqueous solution is in contact with both sides of the with 



«!(*) = -(fcS + l/A^exp^ + l/A^W ^W\RJi(kR)smh(kD/2) 



(k 2 + l/X 2 D y/ 2 e w sinh(fcD/2) + kej cosh(fcD/2) 



2ir\a\kRJ 1 {kR) 

° 2( ' ~ _ (k 2 + l/X 2 D y/ 2 e w sinh(fcD/2) + ke T cosh(fcD/2) ( ' 



I 

and the free energy is where T B is defined from 

2tt 2 , „ „ „,D R , . , 

F = \<r\ 2 R 2 DF B (-, — ) , (A13) 

£j it Ad 

I 



x . 



dkJl{k) 



sinh(te/2) 



cosh(fca;/2) 



„(fc2 +y 2)l/2 



sinh(fcx/2) 



k(k 2 + y 2 y/ 2 (|* + V ^ F) ) sinh(te) + (k2+ 2 y \ )1/2 cosh(fca:) 



(A14) 



In the low screening limit R << Xd, this function be- 
comes of the general form Eq. IC1I with A = — 1) 
and -B = i(| i - + | I — 2), thus the free energy becomes 

F = ^\ a \ 2 R 2 (l ^ ih HR/D) + / a }) , (A15) 
£/ irR 

where the coefficients are given from 



h = — 

£l 

h = ^ln(8/e 1 / 2 ) + (^-l)(^ + ^) 



(A16) 



si 



n. + 1ml _ 2 



2 + £j- + ^ 

Era Si 



<LL_ + e +2 > 

Em El 



£i. _|_ Sm. _ 2 



These expressions are actually identical with Eq. lA7l with 

a = b = ^. 



The high screening limit (Xn << R) corresponds to 
y — > oo. In this limit, T B takes a particularly simple 
form 



T B {x,oo) = 



xy Jo 



dk 



J 1 2 (fc)e / _ lei 

kSqn 



xyey 



leading to the free energy 



2tt 2 

F = 2 —\a\ 2 R 2 X 



D • 



(A17) 



(A18) 



This result is significantly different from the low screening 
limit, as the free energy is independent of the separation 
between the two disks D, which has been replaced by the 
Debye length and the dielectric constant is no longer the 
dielectric constant of the membrane but that of the water. 
This expression is twice the free energy of the diffuse layer 
Eq. [15] and implies that the system minimizes the free 
energy by effectively decoupling the two charged disks 
into two diffuse layers, resulting in a zero electric field 
within the bilayer. This result also follows from noticing 
that the "Guoy-Chapman" free energy Eq. IA18I is lower 
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that the "capacitor" formula Eq. [14] whenever 

D 2A_d 

si e w 



(A19) 



which is always satisfied for a sufficiently small A d ■ There 
is another regime defined by R » Xd and Xd » D 
where the "capacitor" formula is still the leading term, 
but its corrections are 0{D/Xd)- This regime is dis- 
cussed in detail in the next problem. 



3. Two oppositely charged disks in solution 

This situation is similar to the previous case, except 
that now both the solvent and the ions can be found in 
between the two plates. The electric potential follows 
from Eq. QT] with A a D — X b D = A^ . The boundary condi- 
tions are given by the usual Eq. [T51 but only one dielectric 
constant appears (e/ ~ Ea — £ w )- The <Zj(fc) coefficients 
satisfy a 3 (k) = — a2(k) and ai(k) = — a 4 (fc), with 



47T , 



ai(fc) = |fj|i?Ji(A:i?)sinh((fc 2 + l/A^) 1 ^£)/2) 

02(A) = - — \a\RJ 1 (kR)exp(~(k 2 + l/X 2 D ) 1 / 2 D/2) 



The free energy is 



£w it ad 



(A20) 



(A21) 



where T is defined by 



T°{x,y) 



dkJ\{k)- 



exp(-(fc 2 +y 2 y/' 2 x) 



x J k(k 2 + y 2 ) 1 / 2 

(A22) 

The low screening limit defined by Xd >> R (y — > 0) 
has already been considered as it is identical as the one 
considered in Sect. ED ( T c (x,0) = T A (x,Q)), whose 
limit is provided by Eq. IA7I 

In the high screening limit (R >> Xd), we consider 
two situations corresponding to either D >> Xd or 
D << Xd- In both cases the T c function Eq. IA22I needs 
to be evaluated for y >> 1, but with the additional as- 
sumptions that either xy » 1 or xy << 1. For xy >> 1 
it is readily obtained 



F°{x,y) = — (1 - exp(-icy)) 
xy 



and the resulting free energy is 



2tt 2 



F = \a\ 2 R 2 X D (l - cxp (-D/Xd)) 



(A23) 



(A24) 



The first term is the diffuse layer free energy of two fi- 
nite, infinitely thin disks surrounded by a diffuse layer 
of counterions on both sides. In this case, each disk has 



a contribution that is exactly 1/2 of Eq. [T5J The sec- 
ond term represents the Coulomb attraction of the two 
disk-charges, which is screened by the electrolytes. 

The limit xy « 1 is more complex, and it is obtained 
from the following function 



Qj{x,y) = I dkjf(k) 



l-exp(~{k 2 +y 2 ) 1 / 2 x) 



._> M l -exp(-(fc 2 + y 2 ) 1 / 2 x) 



dkJf(k) 



(A25) 

which is related as d(x g ^ C) = 2Qj (see Eq. [£22). The 
first integral is expanded as 



dkJl{k) 



l-exp(-(fc 2 + y 2 ) 1/2 a;) 



(A26) 



£ fjt(k) - xy£ d p 2 (k)(l + k 2 jy 2 fl 2 
K Jy Jo 



The second integral is computed from 



] dk.m l -^-tk 2 + y 2 y> 2 x) 



(A27) 



dk 



(1 - sin(2fc)) exp(-(fc 2 + y 2 ) 1/2 x) 



Ti" Jo k 2 

= — ^(l-sin(2yu))exp(-^(l + M 2 ) 1 / 2 ) 
■ny J x u z 



= —(l + 0(xy)) . 
Try 

where the asymptotic expansion of the Bcsscl functions 
for large values of the argument have been used. Com- 
bining equations Eq. lA26l and Eq. lA27[ and noticing that 
the terms of order 1/y cancel, it follows that 

Qj(x,y) = ^^ + 0((xy) 2 ,(xy)/y) , (A28) 

which leads to the small x-expansion for J- c , 



2.x xy 



x 2 



or 



^ c (x,y) = 1 



xy 



(A29) 



(A30) 



The free energy in the regime defined by D << Xd and 
R » Xd is therefore 

F= ^L\ a \ 2 R 2 D{l ^JL ^ (A31) 



ei 



2Xd ' 



that is, the first correction is independent of the domain 
radius R, consistent with the short-range nature of the 
screened potential. 
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APPENDIX B: INTERACTION OF DOMAINS where Qi are functions defined explicitly 
ACROSS A BILAYER 



In this appendix we consider two uniform dipole den- 
sities separated a distance L, as shown in the middle 
picture of Fig. |4j In this situation we consider that both 
dipole densities are surrounded by a media with the same 
dielectric constant and no salt. The electric potential can 
be obtained from the coefficients Eq. IA1I into Eq. QT] and 
using the superposition principle. Introducing this po- 
tential into the free energy Eq.[5]the following expression 
follows (assuming L>D). 



F = Fq + Fq — 27t(<ti — <T2) 2 -D 2 i21n( 



8R S 



D 2 RH{ 



D 
(Bl) 



where r c 



Fq is the capacitor free energy Eq. [TJ] for a sur- 
face charge Oj and the function Ti is independent of R. 
If the dipole densities are the same o\ = 02, the interac- 
tion free energy cost is linear in the radius R and can be 
incorporated into the line tension. 



Q x (z) = ^((2n + l) 2 ln(2n+l)-4n 2 ln(n))z n 

n=l 

= £)((n+l)lQg(n+l)-nlog(n))z n , (C3) 



n=l 



which arc both convergent for \z\ < 1. It follows from 
Eq. IC2I that if the coefficient in Eq. IC1I satisfies B > 
— 1/2, the argument of Q is \z\ < 1, thus ensuring the 
convergence of the results. 
APPENDIX D: TWO RELEVANT INTEGRALS 



The following results are used throughout the calcu- 
lation and are quoted without further derivation. These 
integrals are calculated to the first non-trivial order in an 
expansion in powers of its argument x. 



APPENDIX C: EXPANSION FOR T 

In many of the situations we consider in this paper, we 
have to find the expansion for small x of a function 



T(x) 



2 J\ (k)(l - exp (-kx)) x 



dk 

K 

1 + A(l - exp {-kx)) 

1 + 5(1 - exp (-2kx)) ' 



(CI) 



where A, B are parameters that depend on the particular 
problem. Upon expanding the denominator in powers of 
exp (— kx) and using the integrals Eq. ID1I and Eq. ID21 
the following expression follows 



x 

~ 2^ 
2 1 



2tt 1 



2—A(l 

7T 



2B) 



B 



B 



-4Ag 2 ( 



B 



l + B 



)) 



x x 
2 ~ 2^ 



f°° dk 

I — 7? J?(fc)(l — exp (— kx)) exp(— nkx) = 
Jo k 

(2n + 1) ln(8/(e 1 / 2 ^)) _ („ + l) 2 \ n (n + 1) + 

+n 2 ln(n)](Dl) 



f°° dk 

/ — 7T Jf(k)(l — exp (— kx)) 2 exp(— nkx) 
Jo k 2 



(l + 4S-2A)ln(8/(e 1/2 x)) (C2) = 



x 
2~k 



2 ln(8/(e 1/2 a;)) - (n + 2) 2 ln(n + 2) + n 2 ln(n) 



(D2) 
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